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ABSTRACT 



We calculate factorized cross sections for lepton pair production mediated by a 
virtual photon in hadron-hadron collisions using the jet calculus scheme, in which a 
kinematical constraint due to parton radiation is taken into account. This method 
guarantees a proper phase space boundary for subtraction terms. Some properties 
of the calculated cross sections are examined. We also discuss matching between the 
hard scattering cross sections and parton showers at the next-to-leading logarithmic 
(NLL) order of quantum chromodynamics (QCD). 

§1. Introduction 

In the evaluation of hadron-hadron scattering cross sections, logarithmic con- 
tributions due to collinear parton production are subtracted from the hard scatter- 
ing cross section, and they are absorbed into the parton distributions of hadrons. 
However, the leading-logarithmic (LL) order of quantum chromodynamics (QCD) 
is insufficient to evaluate these processes. Thus, the next-to-leading logarithmic 
(NLL)-order contributions should also be taken into account. The next-to-leading- 
order (NLO) calculation is also necessary in order to remove theoretical ambiguities 
due to the factorization procedure as well as the choice of the factorization scale. 1 ) 

In actual calculations, Monte Carlo methods are powerful tools for the evaluation 
of exclusive processes. So far, parton shower models for initial-state parton radiation 
have been mainly constructed on the basis of the LL order of QCD, 2 ) although some 
refinements, such as angular ordering conditions for soft gluon radiation, are taken 
into account. 3 ^ 

The scaling violation of parton distributions inside hadrons can be understood 
as the sequential evolution of partons in the initial state. Conventionally, the scaling 
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violation of the parton distributions is calculated by solving the renormalization 
group equations in moments. Then these solutions are numerically inverted to yield 
the momentum fractions of partons. 

Alternatively, one can use parton shower models in order to evaluate the scaling 
violation of the parton distributions. One such algorithm has been proposed in Ref. 
4), in which the parton showers are generated at the LL order of QCD using an 
algorithm consisting of a model based on the evolution of momentum distributions. 
In this model, the scaling violation of the parton distributions is generated using 
only information from the splitting functions of the parton branching vertices and 
input distributions at a given energy. It has been found that the method reproduces 
the scaling violation of the flavor singlet parton distributions up to their normaliza- 
tions. This algorithm has been applied to the Drell-Yan process, 5 ) in which the NLO 
contributions for hard scattering cross sections have been included in a Monte Carlo 
model based on LL-order parton showers. 

In order to allow the application of parton shower models to realistic processes, 
the matching problem should be solved, that is, double counting between hard scat- 
tering and parton showers must be avoided. This problem has been studied 6 ) in 
hadron-hadron collisions with LL-order parton showers. However, this method is 
insufficient, because the subtraction scheme dependence of the collinear singularity 
at the NLO cross section is not canceled by the factorization scheme dependence of 
the distribution functions. 

It may be necessary to combine parton showers at the NLL order for the purpose 
of constructing more accurate Monte Carlo generators. As in the analytic calcula- 
tions, it is desirable to combine NLO cross sections with parton showers at the NLL 
order of QCD. 

Recently, the parton shower model based on the evolution of momentum distribu- 
tions was extended to the NLL order of QCD with the modified minimal subtraction 
(MS) scheme 1 ^ 7 ) as well as that with the jet calculus scheme. 8 ) 

In a previous paper, 9 ) we studied factorization schemes for the hard scattering 
cross sections in hadron-hadron collisions, in which the collinear singularities are 
subtracted using the MS scheme. In the conventional MS scheme, it has been pointed 
out that matching between the hard scattering cross section and the initial-state 
radiation is broken in the calculation of exclusive processes. In order to implement 
factorization schemes that are appropriate for Monte Carlo methods using the parton 
shower models to the accuracy of the NLL order of QCD, we took into account a 
kinematical constraint due to parton radiation when considering the subtraction 
terms, which is called the MS scheme. The momenta of partons in the scattering 
processes are conserved in this method. As an example, we have calculated cross 
sections for the Drell-Yan lepton-pair production mediated by a virtual photon in 
hadron-hadron collisions. We found that the conventional MS scheme gives a double 
logarithmic term, which increases the cross section in the soft gluon limit, whereas the 
MS scheme gives only a single logarithmic term. The infrared behavior is reasonably 
stable in the case of factorization with the MS scheme. However, in the collinear 
region after the collinear singularity is subtracted, negative contributions remain 
for the cross section obtained with both the conventional MS scheme and the MS 
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scheme. Such contributions cannot be ignored at the NLL-order accuracy. Event 
generation employed in the Monte Carlo methods with a negative probability may 
not be appropriate, since the strong cancellation between the negative contributions 
from the hard scattering cross section and the positive contributions from parton 
showers may give unstable results. 

To solve this problem, we examine the jet calculus (JC) scheme studied in Ref. 
8) for the calculation of the Drell-Yan process. In this scheme, an additional phase 
space factor of the hard scattering cross section is also subtracted. As shown in §3, 
most of the collinear contributions to the hard scattering cross section at the NLO 
are subtracted in this scheme. 

At the NLO accuracy, the processes qq — > 7*5, qg — > 7*^ and qg — > j*q con- 
tribute to the hard scattering cross section in actual physical systems. In the follow- 
ing sections, we calculate factorized cross sections for the process qq — ► 7*5 using our 
method. The factorized cross sections for qg — ► j*q are presented in the Appendix 
B. 

In §2, we calculate factorization terms for the collinear singularity by the jet 
calculus scheme. Some properties of the factorized cross sections are presented in 
§3. Section 4 contains a summary and some comments. 

§2. Factorization schemes 

In this section, we calculate factorization terms for the Drell-Yan lepton-pair 
production in quark (g)-antiquark (q) annihilation, 

q(p q )+q(Pq) - 7*(<Z) +g(Pg) - l-(p-) + l + ( P+ )+g( Pg ), (2-1) 

mediated by a photon 7* with the virtuality (p_ + p+) 2 = q 2 = Q 2 , where a gluon 
(g) is radiated in the final state. Here, pi (i = q, q, g) and p± denote the momenta 
of the corresponding particles. 

In the following calculation, the Mandelstam variables are defined by 

S = (p q +Pq) 2 , t=(p q -p g ) 2 , U = (Pq ~ Pgf , (2-2) 

which satisfy s + t + u = Q 2 . 

In order to obtain a finite cross section for the process qq — > 7*3, we subtract 
the collinear contributions due to the branching processes q — > qg and q — > qg from 
the hard scattering cross section. Although these contributions are compensated by 
the parton showers, the remnant of the subtracted cross section at the NLO cannot 
be ignored at the NLL-order accuracy. 

In actual Monte Carlo simulations with parton shower models, the initial-state 
partons are evolved up to a given energy scale M, which corresponds to a fac- 
torization scale for the separation between the initial-state radiation and the hard 
scattering process. The virtualities and longitudinal momentum fractions of the par- 
tons in the initial state are generated according to nonbranching probabilities (form 
factors) and the splitting functions of parton branching processes, respectively. 4 ) 7 ) 8 ) 
The four- momenta of the partons are constructed from these values. 
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The s-channel momentum s of the hard scattering process is constructed by the 
generated momenta p q and p q using Eq. (2-2). Then the virtuality of the photon 
Q 2 at the NLO is generated from the hard scattering cross section. 

The momentum of the quark r q , defined by r q = p q —p g for the branching process 

lipq) ->q(r q )+g(p g ), (2-3) 

is described by in terms of the momentum fraction z q , the virtuality r 2 and the 
transverse momentum r q T, with p q ■ r q T = p q • r q x = 0, as 

Tq = Pq~ Pg = Z q p q + a q p q + T qT , (2-4) 

where a q = r 2 /s for on-shell gluon radiation (p 2 = 0). Here, we set p\= p\ = 0, 
because the relations — p 2 , —p 2 q <C — r 2 are expected in parton shower generation. 
The squared momentum of a photon Q 2 is constructed as 

Q 2 = (r q +Pqf (2-5) 

to 0(a s ) accuracy of QCD. Here, the quantity f = Q 2 /§ is given by 

f = z q + a q . (2-6) 

According to Eq. (2-6), the phase space < — r 2 < M 2 with z q < 1 for the initial- 
state parton radiation covers the region < — t < (1 — f)s for 1 — tm < f < 1 and 
< — t < M 2 for f < 1 — Tm hi the hard scattering cross section. Here, —t = —r 2 for 
p 2 q = p 2 q = 0. Therefore, we consider the two regions 1 — tm < t < 1 and f < 1 — tm 
separately. Here, we define tm = M 2 /s. 

In the region 1 — tm < f < 1, we subtract the collinear contributions from 
within the phase space to the hard scattering process, namely < — t < s(l — f). If 
we ignore the term a q = r 2 /s in Eq. (2-6) and subtract the collinear contribution 
from the hard scattering cross section in the range < — t < M 2 , with f = z q , the 
matching between the hard scattering cross section and the initial state-radiation is 
broken. In this case, M 2 is larger than the kinematical boundary for —t, given by 
3(1 — f) in the region 1 — tm < t < 1. 

As in the previous paper, 9 ) we define, in 4 — 2e dimensions, the subtraction term 
divided by the Born cross section, ctq{Q 2 , ej3> as 

with 

dS qq ] a s 1 
dz q d{-r 2 ) ~ 2vr T(l - e) 



A-irp 2 



P [F] (z 



(2-8) 



The Born cross section ao{Q 2 , e) is given in Ref. 9). 
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where the coefficient of the subtraction term P qq (z q , e) depends on the factorization 
scheme F. In 4 — 2e dimensions, the strong coupling constant is defined by a s y? e for 
the dimensionless coupling a s and a mass parameter fj,. 
We define the integrated contribution as 

H™(e,f), (2-9) 

where / denotes the region of phase space for the hard scattering process being 
considered. Here, sw^W is a limit of the — r q integration. In our previous paper, 9 ) 
we implemented the MS scheme, in which 

for the branching process presented in Eq. (2-3), with 

P£ ) (z q ) = C F \^ L . (2-11) 

Here, Cp = 4/3 is the color factor, and we define 

(fir)) + = f(r) - 6(1 - f) f dyf(y) (2-12) 

for a function /(f) unregulated at f = 1. 

In this paper, we examine the JC scheme, 8 ) in which the coefficient of the sub- 
traction term is given by 

PjfW) = (P q l C] (z q ,e)) + + 5(l-z q ) e -C F , (2-13) 

where 

P\f\z q ,e) = (1 - z q )^P qq (z q ,e) ~ P$(z q ) - eQ q f\z q ) (2-14) 
for e <C 1, with 

Q[f\z q ) = P$(z q ) log(l - z q ) + C F (1 - z q ). (2-15) 

Here, 10 ) 

P qq (z q , e) = P$(z q ) + eP' qq (z q ) (2-16) 

with 

P' qq {z q ) = -C F {l-z q ). (2-17) 

In this scheme, the factor (1 — z q )~ e from the phase space, which gives a negative 
contribution at the NLO, and the 0(e) term of the splitting function P qq (z q ) in 4 — 2e 



/' 
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dS qq 
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dfd(-rf) 2vrr(l-e) 
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dimensions for the branching process 10 ^ are also subtracted from the hard scattering 
cross section. 8 ^ As shown in §3, most of the collinear contributions to the hard 
scattering cross section at the NLO are subtracted in this scheme. In Eq. (2-13), 
the regularization of the infrared singularity is defined so as to conserve the total 
momentum of the initial-state partons. 

First, we consider the region 1 — tm < f < 1 (I = S). For z q = f(F = JC) and 
the subtraction term for the soft gluon radiation is given by 



£[JC](S)( 



M 2 



C F 



<-■-)■■ I d{-r z q ) 
11 



1-f - — J(l-f) 

4 ' 



(2-18) 



Similarly, with z q = f + (— r q )/s (F = JC') and = 1 — f, we obtain 



f[JC](S)/ 
± qq v 



(l-f)S 



Pjf ] (f+(-r 2 )/s,e) 



C F 
3 



+ 



1 1 

~€~Cf~ 
1 



n N /l0g(l -f) 

+ 2(2 + f 2 ) ( 7 . ; 

1 — T 



2(1-*)- 



(2-19) 



The calculation method employed in this region is explained in Ref. 9). 

The factorization term in the region fo < f < 1 — tm (I = C) is given by 



Al 2 



d(-r 



C F 



c^J?^) [~ e + 1o s*m(i - f)J + 1 - f 



(2-20) 



for the JC scheme. Here, we define fo = Qq/s, where Qq is a lower limit for the 
virtuality of the photon to be observed. 
For the JC scheme, we obtain 



= C F 



M 2 







d{-r 2 q ) 



Plf\f+{-rl)/s,e) 



1 ^(*)(- + lo g f^^) + (l + f)lo gr ^- 
V -e 1-t-tmJ 1-t-tm 



lCp~ qq 

+ 1-T-TM 



(2-21) 



*' The function Q^^(z) corresponds to AF\ 3 ^(z) in Ref. 
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The calculation in the collinear gluon region is given in Appendix A. 
By replacing q with q in the above equations, we obtain the subtraction terms 
for the antiquark legs, which are the same as those for the quark legs. 

§3. Factorized cross sections 

First, we calculate the contributions of the soft gluon radiation, which are impor- 
tant for evaluating the NLO contribution for the Born cross section, namely s — Q 2 . 
The factorized cross section to the soft gluon region 1 — tm < t < 1 (I = S) m the 
factorization scheme F is defined by 



dd [F](S) 
uu qq 

df 



— a (Q ,0)K qq " >(t) 



(3-1) 



with 



Ki^ S \f) = K qq (r,e) - F q ^(e,t) - F^ S \e,f) 



(3-2) 



and 



^o(Q 2 ,o) 



4ira 2 
3N C Q 



e 



2 1 



(3-3) 



Here, Nc = 3 and the electric coupling constant of the quark is defined by e q a for 
-e->0. 

The differential cross section for gluon radiation with virtual loop contributions 
is given by 9 ) n 



da q q _a s 2 r(l -e) 



AlTfi 2 



(3-4) 



with 



K qq (r, e) 



C F 



2 1 



log(l ~ r) 
1 -f 



(3-5) 



From Eq. (3-2), we obtain 



qq 



2(1+7* 



2 ^ fl og(l-f) 
1 -f 



2 



Pf q \r)) log r M + 



-2(l-f) 



(3-6) 



*' Our common factor in the hard scattering cross section and the factorization term is different 
from that defined in Ref. If). Therefore, no Pqq(j) logf term appears in Eq. (3-5). 
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with the JC scheme and 



-4 



(1-f) 



+ -- + -7T^U(l-f 



(3-7) 



with the JC' scheme. 

The factorized cross section in the JC scheme depends on the factorization scale 
M. With the JC' scheme, the factorized cross section has no M dependence, because 



the subtraction term is integrated over the range < 



< s(l - f). In this 



region, the virtuality of the photon Q 2 is constructed from the momenta of the 
partons generated by parton showers, in accordance with Eq. (2-6), within the 
region —r 2 < M 2 for z q < 1 . Therefore, the factorization scale for the distribution 
functions in this region is also given by M 2 . 

In order to evaluate the soft gluon contribution, we integrate Eq. (3-1) over 
the range (1 — r] s )s < Q 2 < s, with fixed s, where s is generated by the parton 
showers. Here, r] s is a cutoff parameter satisfying i] s < tm- The energy scale of the 
running coupling constant in this region can be chosen as Q 2 . In order to simplify 
our analysis, we evaluate the soft gluon contribution with the coupling constant 
a s(Q 2 ) — o: s (s). The running coupling constant for the accuracy of the NLL order 
is normalized as a s (M|) = 0.114 at the Z° boson mass. 12 ) 

The integrated cross section is given by 



with 



4f s >(i,i-, fc )^^„( S ,o)4P'(i,i-,,) 



(3-8) 



(3-9) 



The calculated results for the range 1 — r/ s < f < 1 are given by 
4f 1(5) (1, 1 - Vs) = C F [2 {SP-( Vs ,0) - ^(log^ - 2) + log 2 Vs + log(l - n,)} 

-2 I - log(l - r) s ) - r) s + - + 21ogr/ s 1 logf M - ^ + ^vr 2 ] (3-10) 



with the JC scheme, and 
r [JC'](5). 



I q J y J ^(l, l-r ls ) = C F - 4SP. ( Vs , 0) - 2 log 2 r,s + 3 log 



-2 + r 



1 - rj s 



with the JC' scheme. Here, we define 



SP^(a,b)= r df^^. 

Jb 1 - t 



(3-11) 



(3-12) 



9 



We calculate the rj s dependence of the function 



R [F]is \vs 



<T (S,0) 



2tt 



S, /^ 1(5) (l,l 



'qq 



1)s 



(3-13) 



using Eqs. (3-10) and (3-11). 

In Fig. 1, the rj s dependence of the contribution obtained with the JC' scheme is 
plotted by the solid curve. The result obtained with the JC scheme for i] s = tm < 0.5, 
which corresponds to the contribution integrated over the range 1 — tm < f < 1, is 
represented by the dash-dotted curve. The rj s dependences obtained from the two 
schemes become similar, because part of the term log rj s log tm is canceled by the 
log 2 i] s term in the JC scheme, represented by Eq. (3-10). 

The dashed curve represents the r/ s dependence obtained with the JC scheme 
at tm = 0.5. The function I q ^^ near the threshold (rj s <C tm) increases as 
ijjq ^ ~ 2Cf log 2 rjs, whereas the results obtained with the JC scheme for rj s = tm 

and with the JC' scheme decrease as ^ ~ —2Cf log 2 r] s in the soft gluon region. 

The error in the calculation is approximated by 1 < a s ((l — r] s )s)/a s (s) < 1.05 
for < rj s < f M = 0.5 at v 7 ! = 200 GeV. 

As shown in Fig. 1, the soft gluon contribution depends on the kinematical 
boundary for the subtraction term. The result integrated over the region 1 — tm < 
f < 1 with the JC scheme behaves similarly to that with the JC' scheme. Therefore, 
these two schemes give similar results in the soft gluon region. 



CO 

jr 




Fig. 1. The integrated contributions of the cross sections divided by <t (s, 0), which are defined by 
Eq. (3-13) with Vs = 200 GeV. The solid curve represents the result obtained with the JC' 
scheme. With the JC scheme, the dashed curve and dash-dotted curve represent the results 
when tm — 0.5 and tm = r) a , respectively. 



Next, we calculate the remnant in the region tq < f < 1 — tm (I = C), which is 
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given by 



KS ]iC \f) = K£\f,e)-F^\e,f). 



(3-14) 



The collinear contribution of the hard scattering cross section for the real gluon 
radiation integrated over the range < — t < M 2 ( or < — u < M 2 ) is given by 



L 



(r) 



m2 da - 

d(-i) qq 



a. 



dtd{-i) = ^ do{Q2 > e) W 



4-KfI 2 



K$p(t,e) (3-15) 



with 



K$p{r,e) = C F 



c F qq { ' 



— + log : — 

-e \-T -TM 



+ 1 - f - 2f M 



, (3-16) 



where the differential cross section da q q/df/d{—t) is given in Ref. 9). The calcu- 
lation of the remnant in the hard collinear region is presented in Appendix A. We 
obtain 



K^ C \f)=C F 



L-T -TM 



with the JC scheme and 

K^ C \f) = C F 



(1 + f ) log 



(1 



1 - f - f M 



2tm 



TM 



(3-17) 
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with the JC' scheme. 

Here, we calculate the f dependence of the function 



R^ c \t) 



1 



da 



[F](C) 



<m 



<5"o(s,0) df 



a s (Q 2 ) T AF](C) ( ^ 
2vrf qq 1 ' 



(3-19) 



using Eqs. (3-17) and (3-18). 

In Fig. 2, the f dependence of the contribution obtained with the JC' scheme 
is plotted by the solid curve. The result obtained with the JC scheme is represented 
by the dashed curve. Most of the collinear contribution is subtracted from the cross 
section in the JC' scheme, whereas a positive contribution remains near f ~ 1 for 
the result obtained with the JC scheme. 

We also show the calculated results obtained from the MS and MS schemes by 
the dotted curves. These schemes were studied in our previous paper. 9 ) Explicit 
expressions for these schemes are given in Appendix A. 

In Fig. 2, the hard scattering contribution integrated over M 2 < —t < 1(1 — 
f ) — M 2 is represented by the dash-dotted curve. 

As shown in Fig. 2, the remnant of the hard collinear contribution for the 
JC' scheme can be safely neglected even at NLL-order accuracy. We can avoid the 
large cancellation between the hard scattering cross section and the parton showers. 
Therefore, we expect reasonably stable results in Monte Carlo calculations with this 
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Fig. 2. The f dependence of the factorized cross sections with \/i = 10 3 GeV and M — Q — 200 
GeV for the process qq — > 7*g. The solid curve represents the result obtained with the JC' 
scheme. The result obtained with the JC scheme is represented by the dashed curve. The 
calculated results obtained from the MS and MS schemes are shown by the dotted curves. The 
hard scattering contribution integrated over M 2 < — t < s(l — f) — M 2 is represented by the 
dash-dotted curve. 



scheme. On the other hand, the negative contributions of the NLO cross sections 
obtained with the MS and MS schemes should be canceled by the contributions from 
the parton showers. 

Although the parton showers already include the collinear contributions of the 
processes g — > qX and g — > qX, as well as those of q — > qX and q — > qX, at 
NLL order, we have to add the NLO contributions for the processes qg — > j*q and 
qg — > 7*g in the hard scattering region. The calculation of the factorized cross 
sections is rather straightforward. We present the NLO contributions for the process 
qg — > 7*5 in Appendix B. The numerical results of the factorized cross sections for 
this process are plotted in Fig. 3, in which the notation is the same as that in Fig. 
2. In this case, the remnant of the hard collinear contribution with the JC' scheme 
can also be neglected at NLL-order accuracy. 

§4. Summary and comments 

In this paper, we have studied collinear factorization for hard scattering cross 
sections in hadron-hadron collisions, in which the initial-state radiation is generated 
by parton shower models. 

We calculated hard scattering cross sections for the process qq — ► j*g with the 
jet calculus scheme, in which a kinematical constraint for the gluon radiation is taken 
into account. This method is called the JC' scheme. In this method, this kinematical 
constraint guarantees a proper phase space boundary for the subtraction terms. We 



12 



0.15 



i 

o 
E 



0.05 




-0.05 



Fig. 3. The f dependence of the factorized cross sections with \/i = 10 3 GeV and M — Q — 200 
GeV for the process qg — > 7*9 ( gc? — > 7*5 ). The notation is the same as that in Fig. 2. 



also calculated the hard scattering cross section for this process without using the 
kinematical constraint (JC scheme). In the JC scheme, matching between the hard 
scattering cross sections and the initial-state radiation is broken for 1 — tm < f < 1. 

The subtracted cross section integrated over the range 1 — r] s < f < 1 obtained 
with the JC scheme increases as ~ 2Ci?log 2 ?? s in the soft gluon region (r] s <C tm), 
whereas the results obtained with the JC scheme for i] s = tm and with the JC' 
scheme decrease as ~ — ICp log 2 t] s . The kinematical constraint due to the parton 
radiation changes the properties of the NLO cross sections in the soft gluon region. 
If the soft gluon contributions are resummed as 



<t (s,0) + a [ q ^ ]{S \l, 1 - Vs) - <5o(s,0)exp 



2vr 11 \ ' 7,3 > 



, (4-1) 



the cross sections obtained with the JC' scheme and with the JC scheme for rj s = tm 
vanish at the soft gluon limit, which is a desirable property for Monte Carlo methods. 

In the hard collinear region f < f < 1 — tm, after the collinear singularity is 
subtracted, negative contributions remain for f ~ 1 for the cross section obtained 
with both the MS scheme and the MS scheme. Such contributions may not be 
appropriate for event generation employed in the Monte Carlo methods with NLL- 
order accuracy, because large cancellation may occur between the hard scattering 
contributions and those from the initial-state radiation, whereas most of the collinear 
contributions are subtracted from the cross section obtained with the JC' scheme. 
Therefore, with the JC' scheme, we expect to obtain reasonably stable results in the 
numerical calculations. 

Finally we comment on the factorization scheme dependence between the initial 
parton radiation generated by the parton showers at the NLL order and the hard 
scattering cross sections at the NLO. The parton evolution depends on the sub- 
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traction scheme used for the mass singularity, whose appearance is due to collinear 
parton production. 

In parton shower models, infrared unregulated splitting functions are used for 
evaluation of the nonbranching probabilities. For gluon radiation from a quark leg, 
the difference between the splitting function obtained with the MS scheme and that 
obtained with the jet calculus scheme is given by /3oQ qq i (z q )/2 ~ (3oQ qq (f)/2 at 
NLL-order accuracy. x ) Here, (5q = 11 — 2/3Nf with Nf flavour quarks. 

In order to cancel the factorization scheme dependence, we use the relation 

P$ C]il \z) = P^%)-^QW(z) (4-2) 

with 

Qf\z) = P§\') Iog(l -z)- Pljiz) (4-3) 

for i,j = q,q,g at the NLL order in order to obtain nonbranching probabilities. 8 ) 
Here, P^^ 1 \z) and p\^ <yl \z) are the splitting functions calculated with the MS 
scheme 13 ^ and JC scheme, 8 ) respectively. The functions Plj(z) are the O(e) terms of 
the splitting functions in 4— 2e dimensions. 10 ) In this case, the soft gluon contribution 
of the splitting functions may yield large higher-order contributions for z ~ 1. In 
particular, the functions Q^ qq C \z) and Q^gg"\z) behave as ~ log(l — z)/(l — z) for 
z — > 1. Such terms can be absorbed into the strong coupling constant 8 ) as 

2vr 

with Kf- = e^'^ 2 ) / '^ ) ^)(l - z)K 2 ~ (0.4 ~ l)k^, where kr is the transverse mo- 
mentum of the generated parton with spacelike virtuality —K 2 < 0. 

The jet calculus scheme presented in this paper may be useful for the construc- 
tion of a more accurate Monte Carlo algorithm, in which the parton radiation at 
the NLL order as well as the hard scattering cross sections at the NLO of QCD are 
taken into account. 

In future work, we shall construct a Monte Carlo algorithm with the JC' scheme 
using the method presented in this paper. 
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Appendix A 

In this appendix, we derive the remnant in the hard collinear region (/ = C) for 
the process qq — > j*g . 

The factorization terms are given by 



^ ](C) (^)=(i- + logf M )p[f](f, e ) 



for F = MS,JC, where 



and 



Pf q C \r,e) = P^{r)-eQ^{r). 



qq 



The functions Pqq and are given by Eqs. (2-11) and (2-15) in the main text, 

respectively. 



For MS andJC', we obtain 

F^ C \e,r) = C F \±m(r) ( ±- + log 



C F 



(1 - f)f M 

l-T-TM 



and 



1 - f 

+ (1 + f)log- f M 

L-T-TM 

Flf^\e,f) = F^\e,f) + Q[fHf), 



respectively. 

The remnant in the hard collinear region is defined by 

(C) 

where if^ ; (f, e) is given by Eq. (3-16) in the main text. We obtain 

Kft m (r) = C F \-LpV\f ) log -^t- + (1 - f - 2f M ) 
Of 1 — r — tm 



with the MS scheme and 

K mic) {f) = K mio) if) + CF 



2 1 - f - f M 

log — ; h r M 



1 - f 



1 -f 



with the MS' scheme. The results obtained with the JC scheme and JC' scheme are 
respectively given by 
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and 



which are presented in Eqs. (3-17) and (3-18) in the main text. 

The contribution of the hard gluon radiation (7 = H) integrated over the range 
M 2 < -i < 3(1 - f) - M 2 is given by 



K ( q H q \f) = - d{-t)C F 
s Jm 2 



2C F 



^4 0) (r)log 



1 ~ T - f M 

TM 



(1 



2f M ) 



for f < 1 — 2tm- 



Appendix B 



In this appendix, we calculate the cross section of the Drell-Yan lepton-pair 
production in antiquark (g)-gluon (g) scattering, 

q(pq) + g(p g ) 7* («) + Q(p'q) r(p_) + i + ( P+ ) + q(p' 9 ), 

mediated by a photon 7* with virtuality (p- + p + ) 2 = q 2 = Q 2 . Here, pi (i = q,g) 
and p± denote the momenta of the corresponding particles. 

The momentum of the quark r q , defined by r q = p g —p q for the branching process 
dip 9) ~^ li r q) + i s described by using the momentum fraction z q , the virtuality 

r 2 and the transverse momentum r q x, with p q ■ r q T = p g • r q x = 0, as 

r q = z q p g + a q p q + r qT , 

where a q = r 2 /s for p'? = with p 2 = 0. 

In the MS scheme, the coefficient of the subtraction term is given by 



p[MS] ( 



Z q^) 



Pqgi Z q) 



with 



P^{z q )=T R [z 2 q + (l-z q f 



Here, Tr = 1/2 is the color factor. 

In the jet calculus scheme, we define 8 ) 

Pqg \ z qi e ) = Pqg i z q) 



eQq>g C] iZq) 



with 



Q [ qg C] {Z q ) = P$ iZq) log(l - Z q ) + Z q (l - Z q ) . 
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For F = MS,JC, the subtraction terms for 1 — tm < f < 1 (I = S) and 
fo < f < 1 — tm {I = C) are given by 

For F = MS', JC', we obtain 



Fjf ](^)( e ,f) = r fl [-L^pW (f) + i-p(0)(f)log(l - f) 
+(l-f)(3f -1) 

and 

Flf^ s \e,r) = Flf^\e,f) + Ql^ (f) 
for 1 - f M < f < 1 (J = S), and 



FlfW) {e ,f)=T R 



1 1 



,^p(°)(f) + — Pg)(f) log f M 



and 



+2(2f - l)f M + t 2 m 



F^ c \e,f)=F^\e,f) + Q^(f) 
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for the region fo < f < 1 — tm (/ = C), respectively. 
We obtain the subtracted contributions 



< 1(J) (r)=^ ) (f, e )-^« (e ,f) 



with 



1 f( 1 ~' ? )s 



(-*)(-«) 



for I = S, and 



Kf g \f,e) = -l d(-t) 
s Jo 



Tr 



1 



W - + log(l " f ) 2 - 1 + £ + f - £f 



— e 



M 2 



(-f)(-u) 



= T R 



1 

Tr 



P§\t) f — + logf M (l - f) ) + 2f (1 - f + T M ) + ;jf 



'Af 



for I = C. The differential cross section for the process qg — ► 57* is given by 11 ) 



dfd(-t) 2vrs 



(-*)(-&) 

47TS/X 2 



with 



K^(f,-i,e)=T R 



where 



with 



P qg (t,e)=P$(t)+eP qg {t) 



P qg (r) = -f(l-r). 

Here, the Born cross section ao(Q 2 ,e) is given in Ref. 9). 
For I = S, we obtain 



K mw {f) = TR 



1 



(i-f) 



r\2 



,^pW(f)(log^— — 1 ] + 



3 
2 



2 



^ 0) (f)(log(l-f)-l) + ^-3f + V 



with the MS scheme and 

K mw {f) = TR 

with the MS' scheme. 
For I = C, we have 



K [MS](C) (f } = Tr _p(0) (f ) lQg(1 _ f ) + 2f (1 _ f + f M ) + _^ 



and 



K [ ™ S ' m (t) = T R \^r R PS\r) log(l - f) + 2f(l - f - f M ) + ^r M (4 - f M ) 



The calculated results with the JC scheme and JC' scheme are respectively g 



by 



< i( %)=4r ( %)-« ] (T) 



and 



<' 1(J) (f)=< S1 «(f)-Q[f(f), 



for / = S, C. The numerical results are presented in the main text 

) 

( 1 f) 2 



For f — ► 1, the functions (t) behave as 



^s ]( 5) (f) ^ TRlog !i^ idf ]( %) _ Ti? iog(i -f) 
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for F = MS, MS , and 

^cK,) (f) ^ rRlog (L z i) ) 4f ](5) (f)-o 

TM 



for F = JC, JC'. The integrated results for the range 1 — rj a < f < 1, 



l—ris 



T 



vanish for r/ s — > 0. 

The contribution of the hard gluon radiation (I = H) integrated over the range 
M 2 < — i < s(l — f ) is given by 

1 /■«(!-*) 



s JM 2 



r, 



^J? (r) log ^-^ + 2f (1 - f - f M ) + J{(1 - f ) 2 - fU 



for f < 1 — TM- 

The subtracted cross sections for the quark legs are the same as those for the 
antiquark legs. 
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